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ABSTRACT 

We present the results of X-ray spatial and spectral studies of the "mixed-morphology" supernova 
remnant IC 443 using ASCA. IC 443 has a center-filled image in X-ray band, contrasting with the shell- 
like appearance in radio and optical bands. The overall X-ray emission is thermal, not from a synchrotron 
nebula. ASCA observed IC 443 three times, covering the whole remnant. From the image analysis, we 
found that the softness-ratio map reveals a shell-like structure. At the same time, its spectra require 
two (1.0 keV and 0.2 keV) plasma components; the emission of the 0.2 keV plasma is stronger in the 
region near the shell than the center. These results can be explained by a simple model that IC 443 has 
a hot (1.0 keV) interior surrounded by a cool (0.2 keV) outer shell. From the emission measures, we infer 
that the 0.2 keV plasma is denser than the 1.0 keV plasma, suggesting pressure equilibrium between the 
two. In addition, we found that the ionization temperature of sulfur, obtained from H-like Kq to He-like 
Ka intensity ratio, is 1.5 keV, significantly higher than the gas temperature of 1.0 keV suggested from 
the continuum spectrum. The same can be concluded for silicon. Neither an additional, hotter plasma 
component nor a multi-temperature plasma successfully accounts for this ratio, and we conclude that 
the 1.0 keV plasma is overionized. This is the first time that overionized gas has been detected in a 
SNR. For the gas to become overionized in the absence of a photoionizing flux, it must cool faster than 
the ions recombine. Thermal conduction from the 1.0 keV plasma to the 0.2 keV one could cause the 
1.0 keV plasma to become overionized, which is plausible within an old (3x10 yr) SNR. 

Subject headings: conduction — plasmas — radiation mechanisms: thermal — supernova remnant — 
supernovae: individual (IC 443) 



1. INTRODUCTION 

Supernova remnants (SNRs) have been traditionally 
classified into three categories: shell-like, Crab-like (Pleri- 
onic; its emission mainly originates from a pulsar nebula), 
and composite (shell-like containing a Plerion). These 
SNRs usually have similar X-ray and radio morphologies. 
In addition to these, there is a recently established class 
characterized by a center-filled morphology and thermal 
emission in the X-ray band, contrasting with a shell-like 
structure in the radio and optical bands. Rho & Pe- 
tre (1998) termed these "mixed-morphology" SNRs (MM 
SNRs) . The distinctive characteristics of this category are: 
(i) the dominant X-ray emission mechanism is thermal de- 
spite an X-ray morphology similar to Crab-like SNRs; (ii) 
the emission arises primarily from the swept-up interstellar 
material (ISM), not from ejecta; and (iii) they are usually 
near or interacting with molecular clouds or other large 
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ISM structures. 

Rho & Petre (1998) identified two formation mecha- 
nisms that might explain the properties of MM SNRs. One 
is the cloud evaporation model, in which the interior X- 
ray emission arises from the gas evaporated from shocked 
clouds (White & Long 1991). For this model to match 
the properties of MM SNRs, the clouds must be numer- 
ous, small in size, and of sufficient density to survive in 
the passage of a strong shock. They pass into the hot in- 
terior of the remnant, gradually evaporating and emitting 
X-rays. The other is the thermal conduction model (Cox 
et al. 1999; Shelton et al. 1999). In this model, the shell 
of an expanding SNR cools below ~10 6 K and its X-ray 
emission becomes too soft to be detected through the ab- 
sorbing ISM. The thermal energy stored in the shocked 
interior gas is conducted outward, producing the observed 
uniform interior temperatures. This model requires suffi- 
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cient interior density for conduction to be viable. 

IC 443 (G189. 1+3.0) is one of the prototypical MM 
SNRs. It has a ^45' diameter and is located at the dis- 
tance of 1.5 kpc (Fesen 1984). It is strongly interacting 
with a foreground molecular cloud and an H I cloud to the 
east (Cornett, Chin, & Knapp 1977; Giovanelli & Haynes 
1979). It is a candidate for the CGRO EGRET 7 -ray 
source 3EG J0617+2238 (Hartman et al. 1999). Its X-ray 
morphology was observed by Einstein to be center-filled, 
which is highly atypical of a SNR in the adiabatic phase, 
with little correlation with the radio or optical image. The 
spectrum of the brightest area was well described by ei- 
ther a 1.5 keV non-equilibrium ionization (NEI) plasma 
or 0.2 keV and 0.95 keV ionization equilibrium plasma 
components (Petre et al. 1988). Asaoka & Aschenbach 
(1994) analyzed ROSAT all-sky survey and pointing data, 
and found another partially overlapping SNR of ~T0 5 yr 
age, which they called G189.6+3.3. They reproduced the 
spectra for the regions excluding the part overlapping with 
G189.6+3.3 by a two-temperature (1 keV and ^0.3 keV) 
ionization equilibrium plasma model. 

Using Ginga, Wang et al. (1992) found a hard X-ray 
component in the spectrum of IC 443 extending up to 
20 keV. They suggested that this component originates 
from a very high temperature (>10 8 K) plasma and thus 
the remnant's age is about 10 3 yr instead of the previous 
estimate of 10 4 yr. However, ASCA and BeppoSAX re- 
vealed that this hard component is non-thermal and orig- 
inates in two small regions along the southern edge of the 
radio shell where the SNR-molecular cloud interaction is 
strongest (Keohane et al. 1997; Bocchino & Bykov 2000). 
Using Chandra data, Olbert et al. (2001) determined that 
one of the regions is a synchrotron nebula powered by a 
compact source (probably a neutron star), whose kick ve- 
locity is consistent with an age of 3xl0 4 yr. 

In this paper, we first present the ASCA GIS images 
and softness ratio map (§ 3.1). This reveals the existence 
of a soft, shell-like structure. We compare the SIS spectra 
of the north region near the shell with that of the central 
region (§ 3.2) and then perform the spectral fitting (§ 3.3). 
Finally, we discuss a possible mechanism that can account 
for our results on the ionization state of IC 443 (§ 4). 

2. OBSERVATIONS 

ASCA has two Solid-state Imaging Spectrometers (SIS0, 
SIS1) and Gas Imaging Spectrometers (GIS2, GIS3) at the 
focal planes of the X-ray telescopes (XRT). For IC 443, 
three observations were carried out. Figure 1 shows the 
fields of view (FOVs) of the observations in the current 
analysis (20' radius region in GIS and entire SIS FOV) 
overlaid on the Palomar Observatory Sky Survey red im- 
age. The relevant information about these observations is 
given in Table 1. We extracted these data from the ASCA 
public archive. We used all available GIS data. During 
the AO-1 phase, the GIS3 pulse height analyzer suffered 
from an on-board processing problem: the least significant 
PHA bits were stuck, resulting in a loss of spectral resolu- 
tion. Therefore the GIS3 AO-1 data used here have only 
8 binned PHA channels. For the SIS, we used only the 
PV phase observation data. There are enough photons 
from this observation to allow investigation of the spatial 
variation of the X-ray spectra with good energy resolution. 



For the most part, the NASA/GSFC revision 2 standard 
data processing criteria were used for screening all the GIS 
and SIS data. The lone exception was for the earth eleva- 
tion angle for the SIS: in order to avoid contamination by 
stray light, we used only SIS0 data taken at angles >25° 
and SIS1 data taken at angles >20°. Since the PV phase 
SIS data show no radiation damage effects, we did not 
apply the Residual Dark Distribution (RDD) correction. 

As IC 443 is located in the Galactic anti-center region 
(l,b) = (189.1°, 3.0°) where the Galactic Ridge emission is 
negligible, we used blank-sky data that include the cosmic 
X-ray background (CXB) and the non- X-ray Background 
(NXB) for background subtraction. The PV phase GIS 
data, with no spread discriminator application, contain 
more NXB events than those in subsequent observations. 
For these data, we used the background datasets that are 
applicable to observations performed when the GIS spread 
discriminator was turned off. 

In the present analysis, we used DISPLAY45 version 
1.90 for the image analysis, and XSPEC version 11.0 for 
the spectral analysis. 

3. DATA ANALYSIS AND RESULTS 

3.1. Image Analysis 

We generated a complete image of IC 443 using GIS 
data. The three GIS2 and GIS3 images were overlaid, 
and an exposure-corrected blank-sky image was subtracted 
from it. We multiplied the simulated XRT efficiency map 
by the GIS grid map to form the vignetting image. Each 
map is smoothed by the corresponding point spread func- 
tion (PSF). We constructed the exposure-corrected vi- 
gnetting image, and divided it into the source image to 
create a vignetting- and exposure-corrected composite im- 
age. Figure 2(a) shows the resultant GIS image in the 
0.7-10.0 keV band with ASCA FOVs. 

In every energy band, the composite image shows IC 443 
to have a centrally-filled morphology. The remnant is 
larger in images below 1.5 keV than in higher energy 
bands. Thus a softness-ratio (F . 7-1.5 kov/^i.s-io.o kev) 
map reveals a shell-like structure, contrasting with the 
center-filled appearance in the X-ray maps (see Fig. 2(b)). 
As shown in Figure 2(c), the soft shell corresponds well 
with the bright optical emission in the northeast and 
southwest. The contours are limited at the FOV edge 
along the northwest boundary, indicating that it extends 
further. The absence of the softness-ratio shell in the 
southeast and its low contrast in the northwest may be 
due to the absorption by the molecular cloud in front of 
IC 443 (Cornett, Chin, & Knapp 1977). 

In Figure 3 we show background-subtracted and 
exposure- and vignetting-corrected SIS images in the 0.5- 
1.0 keV and 1.0-2.0 keV bands. The SIS FOV includes 
both the X-ray brightest region of IC 443 and a part of 
the northeast shell. The peak surface brightness of the 
soft (0.5-1.0 keV) component appears to the north of that 
of the hard (1.0-2.0 keV) component (i.e. closer to the 
rim). 

Therefore, both the GIS and SIS images suggest that a 
soft (<1.0 keV) X-ray component dominates near the shell 
while a hard (>1.5 keV) one dominates in the interior. 

3.2. Difference between Two Regions 
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We therefore searched for a spectral variation from the 
interior to the rim. We used only SIS data for the spectral 
analysis because the SIS has superior energy resolution to 
the GIS (FWHM of 5% at 1.5 keV during the observa- 
tion (Tanaka, Inoue, & Holt 1994)) and a larger effective 
area below 1 keV. We extracted spectra from the two re- 
gions indicated by boxes in Figure 3: a region near the 
northern rim (hereafter North) where the soft X-ray emis- 
sion is stronger, and an interior region (hereafter Center) 
where the hard emission is stronger. Neither region over- 
laps G189.6+3.3, and stray light from the outside of the 
SIS FOV is negligible as shown in Figure 2(a). The spectra 
and their ratio are shown in Figure 4. 

Above 1.4 keV, the shapes of the spectra appear similar 
to each other (see Fig. 4(b)), and are in good agreement 
with the spectrum from a thin thermal plasma with a tem- 
perature of about 1 keV. Below 1.4 keV, on the other hand, 
the North spectrum has relatively stronger emission than 
the Center, and shows the Ne IX emission line at 0.92 keV. 
The intensity ratio of the Center spectrum to the North 
below 1.4 keV is proportional to the energy, as shown in 
Figure 4(b). This ratio would be an exponential function 
of energy if the difference were due to higher column den- 
sity. Therefore, it is difficult to explain this difference as 
a result of column density variation. 

3.3. Spectral Analysis 

In order to characterize the spectral properties and eval- 
uate the difference between the two regions quantitatively, 
we performed model fittings; first with a one-component 
plasma model, and then with a two-component plasma 
model. Since the SIS quantum efficiency changes greatly 
at the energy of the oxygen edge (0.53 keV), we did not 
use the data below 0.6 keV. 

3.3.1. One- Component Plasma Model 

We first applied a simple single-temperature, singlc- 
ionization-timescale, non-equilibrium ionization (NEI) 
plasma model (NEI in XSPEC, Borkowski (2000)) that was 
favored in the previous IC 443 analysis (Petre et al. 1988). 
However, the minimum values of reduced x 2 were 924/118 
for the North and 1337/128 for the Center, even with vari- 
able abundances, which are hardly acceptable. 

We therefore tried an ad-hoc model consisting of ther- 
mal bremsstrahlung, plus 12 Gaussian profiles to represent 
the most prominent emission lines above 0.9 keV. The line 
widths and strengths were all left free. Such models have 
been used previously to model ASCA spectra for other 
remnants (e.g. Miyata et al. (1994); Holt et al. (1994)). 
Both H-like and He-like Ka lines of Mg, Si, and S were 
well fitted with negligible line widths compared with the 
energy resolution of the SIS. 

The flux of each Gaussian profile is proportional to the 
cmissivity of the line component. The line intensity ratio 
of H-like Ka to He- like Ka (hereafter -Rh/Hc) provides a 
measure of the degree of ionization of a given element, and 
depends on both the electron temperature and the ioniza- 
tion timescale. We evaluated the i?H/Ho f° r Mg, Si, and 
S through comparison with the theoretical values (Mewe, 
Gronenschild, & van den Oord 1985) as shown in Figure 5, 
and found that it is impossible to explain the degree of 
ionization of all elements with the same ionization tem- 



perature in each region. This suggests that more than one 
plasma component is required. 

3.3.2. Two-Component Plasma Model 

We next applied a two-component plasma model. The 
continuum components above 1.4 keV of the two spectra 
have a shape consistent with a temperature of 1 keV, as 
mentioned previously. In addition, Si and S appear to be 
ionized enough to apply a collisional ionization equilibrium 
(CIE) model at this temperature (see Fig. 5). We thus 
used the CIE plasma model (VRAYMDND in XSPEC, Ray- 
mond & Smith (1977)) as the high temperature plasma. 
On the other hand, we used the generalized NEI model 
with varying temperature (VGNEI in XSPEC, Borkowski 
(2000)) as the low temperature plasma. All abundances 
were treated as variables, and those of two components 
linked together. This model failed to reproduce the Ka 
lines of Ne X, Si XIV, and S XVI (and K/3 line of Ne X 
in Center), so we added several narrow Gaussian compo- 
nents. 

We were able to fit reasonably both North and Cen- 
ter spectra (except for the energies around the emission 
lines of Si and S) with a model consisting of ~0.2 keV 
NEI and ~1.0 keV CIE plasma with a column density of 
7.4 x 10 21 cm~ 2 as shown in Figure 6. The best-fit parame- 
ters of the model are listed in Table 2. All abundances are 
smaller than the solar values (Anders & Grevesse 1989), 
which suggests that the X-rays originate from the shocked 
ISM. Although these fits are not formally acceptable, the 
result implies that the spectrum can be generally described 
by two plasma components: Figure 6 suggests that the 
spectra below 1.4 keV are well reproduced by the low tem- 
perature (~0.2 keV) plasma, whereas those above 1.4 keV 
show the features of the 1.0 keV plasma. 

Almost all fitted parameters for the two spectra are con- 
sistent with each other within the 90 % confidence range. 
This suggests that the plasma properties in the two regions 
are similar. The sole significant difference is the ratio of 
the emission measures of the 0.2 keV plasma to the 1.0 keV 
one: it is substantially higher in the North (185+gg 2 ) than 

in the Center (48+24)- Hence, we suggest that the differ- 
ence of the intensity ratio between two plasma components 
produces the apparent difference of the spectra shown in 
Figure 4. 

While the two component model accounts reasonably 
well for the gross spectral features, it leaves an important 
detail unresolved. The H-like lines of Ne, Si and S are so 
strong that it is necessary to add Gaussian components to 
represent them. 

In order to investigate the anomalous S line ratio, we fit 
the 2.2-6.0 keV spectra with a model consisting of a CIE 
plasma and three narrow Gaussian components (S XV Ka, 
S XVI Ka, S XV K/3). (We can ignore the contribution of 
the 0.2 keV plasma component at the energies around the 
S lines.) The fits are acceptable as shown in Figure 7. The 
i?H/Hc for S requires an ionization temperature (hereafter 
T z ) of about 1.5 keV in each region, higher than the con- 
tinuum temperature (hereafter T e ) of 1.0 keV as shown in 
Figure 8. The same can be concluded for Si. Ne X emits 
mostly around kT c = 0.5 keV, which might indicate that 
the strong Ne X line comes from an intermediate temper- 
ature (0.2< kT c <1.0 keV) plasma. 
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To account for the higher ionization temperatures of Si 
and S, we replaced the Gaussians with a second, hotter 
(2.0 keV) CIE component, whose temperature is selected 
to maximize the emissivity of the H-like line of S. How- 
ever, this component requires an extraordinarily large S 
abundance; more than 10 3 solar. Such a large abundance 
is not observed even in ejecta, and thus this model is phys- 
ically implausible. Moreover, narrow band images at the 
energies corresponding to Si and S H-like Ka lines show 
no conspicuous (e.g. clumpy) structure and are similar to 
those of other bands (e.g. images for Si and S He-like Ka 
and continuum bands), indicating that the H-like lines are 
emitted from the same plasma that emits He-like lines. 

Another possible source of the anomalous ratio between 
T z and T c is a temperature distribution in the plasma, 
as observed in magnetic cataclysmic variables (Ezuka & 
Ishida 1999). However, as in SNRs such plasma shows 
T z < T c , and thus does not explain the higher T z of Si 
and S. Therefore, we suggest that the Ka lines of H-like Si 
and S are emitted from an "overionized" 1.0 keV plasma 
(i.e., a plasma in which the degree of ionization is larger 
than that expected from the electron temperature assum- 
ing collisional ionization equilibrium). 

4. DISCUSSION 

4.1. The Plasma Structure of IC U3 

ASCA observations have revealed the following new X- 
ray features of IC 443; (1) the X-ray softness-ratio map 
shows a shell-like structure that correlates with the op- 
tical filaments, implying that the hard X-ray emission is 
more centrally concentrated than the soft; (2) the spec- 
tra include two plasma components with temperatures of 
0.2 keV and 1.0 keV, and the intensity ratio of the former 
to the latter is larger in the outer region than that in the 
inner region. 

These results suggest that IC 443 has the plasma struc- 
ture that can be modeled by a central hot (1.0 keV) re- 
gion, surrounded by a cooler (0.2 keV) shell. The mean 
electron densities within these regions are estimated from 
the observed emission measures. The emission measure 
of the 1.0 keV plasma, which is estimated from the CIS 
spectrum of the brightest region within 16' radius (north- 
east part of IC 443 where most X-rays are emitted), 
is ~ 1.1 x 10 58 (d/1.5 kpc) 2 cm~ 3 (d is the distance to 
IC 443). In contrast, that of the 0.2 keV plasma is 
~ 5.4 x 10 59 (d/1.5 kpc) 2 cm~ 3 , more than an order of 
magnitude larger. To estimate the densities of the two 
components, we assume that the 0.2 keV plasma extends 
to the north-east optical shell, and thus has an outer ra- 
dius of 02 — 16'. The radius of the region containing the 
1.0 keV component is assumed to be Q\ — 10', considering 
that the surface brightness of the 1.0 keV at this radius 
is 50 percent of its maximum (from the 2.0-3.5 keV band 
image). The densities of the 1.0 keV plasma (ni) and the 
0.2 keV plasma (712) are 

— °(w)" /2 (^b)" 2 cm ~< (1) 

where a — 9 2 /0i- The cooler exterior shell is denser than 



the interior; the two components appear to be near pres- 
sure equilibrium, and consistent with a typical SNR pre- 
dicted by the Sedov- Taylor solution, despite the fact that 
the X-ray image appears center-filled. 

This structure is also consistent with the prediction of 
the conduction model for mixed morphology SNRs, de- 
veloped for W44 (Cox et al. 1999; Shclton et al. 1999). 
These authors suggest that if a remnant, evolving in a 
fairly smooth ambient medium of moderate density with 
thermal conduction active in its interior, reaches the phase 
that radiative cooling is important, then it shows inner- 
hot and outer-cool regions of X-ray emission interior to 
the cold dense shell. IC 443 has similar characteristics to 
W44, including a partial H I shell, suggesting it is in the 
same evolutionary state. 

4.2. Formation Mechanisms of the Overionized Plasma 

Additionally, we found that the Ka lines of H-like Si 
and S are unusually strong, indicating that the plasma in 
IC 443 is overionized. This is the first detection of evi- 
dence for an overionized plasma in a SNR. In SNRs, we 
are accustomed to consider the plasma undcrionized, as 
in a collisionally-ionizcd, shock-heated plasma, the ioniza- 
tion timescale is typically large compared with the time 
since the bulk of the gas was shock heated. Neverthe- 
less, conditions can (and apparently do) exist in which a 
shock-heated, low-density plasma can become overionized. 
In this section, we discuss possible formation mechanisms 
for such a plasma. 

We can consider two possible ways to produce an overi- 
onized plasma. One is that photoionization causes heavy 
elements to be more highly ionized than their collisional 
equilibrium level. A low-density plasma does not pro- 
duce sufficient radiation to photoionize the heavy ele- 
ments: the plasma needs a bright external radiation 
source. CXOU J061705.3+222127, with the 1-5 keV flux 
is 2 x 10~ 13 erg cnT 2 s" 1 (Olbert et al. 2001; Bocchino & 
Bykov 2001), is too faint to photoionize the 1.0 keV plasma 
to emit the "extra" H-like S flux of ~ 10~ 13 erg cm~ 2 s _1 
shown in Table 2, even if it was at the center of the rem- 
nant in the past. Futhermore, there is no other evidence of 
such a strong source in IC 443 and it is implausible to con- 
sider that this plasma region is illuminated continuously 
by a hidden intense source. 

The other way is that the gas cools faster than the ions 
can recombine. The three possible cooling mechanisms in 
the gas are radiation, expansion, and thermal conduction. 
Shclton (1998, 1999) has shown that inclusion of all three 
processes in hydrodynamic modeling of SNRs expanding 
in a very low density ISM leads to the development of an 
"overionized" plasma in the interior once a remnant has 
reached the radiative stage. The H I shell covering the 
eastern part of the remnant (Giovanelli & Haynes 1979) 
suggests that at least part of IC 443 is in the radiative 
stage. Below we compare the cooling timescales of radi- 
ation, expansion, and conduction with the recombination 
timescale of heavy elements in the 1.0 keV plasma, and 
show that these are consistent with the formation of an 
overionized plasma in IC 443. 

The recombination timescale of an element of atomic 
number Z, which is the characteristic timescale to reach 
ionization equilibrium, is calculated by Masai (1994). In 



ASCA Observations of SNR IC 443 



5 



the case of sulfur at a temperature of 1.0 keV, this is esti- 
mated to be 

z 

Womb « & + a,)" 1 - 9 x 10 11 ' " " 

z=0 



1 cm 3 



-l 



(3) 



where S z and a z represent the rate coefficients for ion- 
ization and recombination from an ion of charge z to 
charge z+1 and z— 1, respectively. This is consistent 
with the expectation that collisional ionization equilib- 
rium is approached for the ionization timescale value n c t ~ 



We estimate the cooling timescale of the 1.0 keV plasma 
via radiation to be 

3 / 2 / A \ V2 

S, 



icooi = — j ~ 4 x 10 i4 . 

-^rad V 10 



1.5 kpc 



(4) 

where V and L rac j are the volume and the luminosity of 
the 1.0 keV plasma respectively. The cooling timescale of 
equation (4) is more than two orders of magnitude larger 
than the recombination timescale of equation (3). 

Rapid adiabatic cooling and the resultant ovcrioniza- 
tion have been discussed by Itoh & Masai (1989) for a 
young remnant breaking out of the circumstellar matter 
into a rarefied medium. However, this is not the case for 
IC 443, which is presumed to be in its Sedov phase (ex- 
pansion of the shocked ejecta is negligible compared to the 
blast wave). The shock temperature decreases with blast 
wave expansion. Even if the effect is taken into account, 
it is found that the electron temperature varies as slowly 
as T e oc t~ 2 / 25 with Coulomb collisions in the post-shock 
region (Itoh 1978; Masai 1994). Therefore, adiabatic ex- 
pansion is unlikely to work, since the timescale is much 
longer than the recombination timescale. 

We next calculate the energy transport rate by ther- 
mal conduction from the 1.0 keV interior to the 0.2 keV 
outer region. We assume a temperature gradient scale 
length, l T = (gradlnT) -1 , of (l-2)xl0 19 cm, the dis- 
tance from the approximate boundary between the 1.0 keV 
and the 0.2 keV components to the edge of the remnant. 
This distance is much longer than the mean free path of 
electrons A c « 10 18 (T c /10 7 K) 2 (n e /l cm" 3 )" 1 cm. The 
total flow of heat to the outer plasma is estimated to 
be 2 x 10 38 (? T /10 19 cm)- 1 (6 1 /10') 2 (d/1.5 kpc) 2 erg s" 1 . 
The energy loss via radiation in the 0.2 keV plasma is 
7xl0 37 erg s _1 from the observed emission measure, com- 
parable to the heat flow. In this circumstance, thermal 
conduction equilibrates the temperatures of the two com- 
ponents on the conduction timescale, 
-l 



t, 



cond 



dlnTr, 



dt 



A; 



2 x 10 



U cm" 3 J VlO 19 cm/ V 



kTx 
1.0 keV 



-5/2 



negligibly small in mature SNRs. We conclude from this 
argument that it is entirely possible in IC 443 and similar 
SNRs for the interior gas to be an overionized, recombining 
plasma. 

Assuming both T z and T are e-folding with their char- 
acteristic timescales, we can estimate the age of the 
1.0 keV plasma (t age ) approximately by the following equa- 
tion, 



Tz ^ exp [ ^age/^recomb] 
T e exp [ ^agc/^cond] 



1.5, 



(6) 



In 



32 



where In A is the Coulomb logarithm. The gas in the 
1.0 keV plasma cools on the timescale given by equa- 
tion (5). This timescale is smaller than the recombina- 
tion timescale, and is comparable to it even with It = 
2 x 10 19 cm. This indicates that the strong gas cooling via 
conduction makes T e lower than T z . While the presence of 
a magnetic field will in principle reduce the heat flux, we 
ignore its effect, since the magnetic field is expected to be 



to be (0.3-1) x 10 4 yr. This age is younger than that esti- 
mated by Chevalier (1999). It indicates at least that the 
1.0 keV plasma is younger than i rcC omb ~ 3 x 10 4 yr; if this 
plasma is older than the age represented by the recombi- 
nation timescale, we would expect T z to be comparable to 

Te- 

Evidence of overionization has not been detected in any 
other SNR. We consider our detection in IC 443 to be plau- 
sible for the following reasons: (1) The interior of IC 443 
was shocked long ago and has had no subsequent source 
of heating or ionization, and thus should be much nearer 
ionization equilibrium than newly shocked gas near the 
rim or in young SNR's; (2) IC 443 has an interior plasma 
hot enough to give rise to strong thermal conduction; (3) 
ASCA SIS is the first detector with sufficient energy res- 
olution to cleanly resolve the H-like and He-like lines of S 
from each other; (4) The signal-to-noisc in the spectra is 
sufficiently high ((4-5) xlO 4 counts in each region) for T z 
and T e to be inferred with small errors. It should be noted, 
however, that no evidence of overionization is detected in 
W44, despite it having similar properties and comparable 
quality of ASCA spectra. The difference between the two 
remnants merits further investigation. 

5. CONCLUSION 

We have discovered a shell-like structure in the softness- 
ratio map and two distinct (0.2 keV and 1.0 keV) plasma 
components in IC 443. Using these results and the emis- 
sion measures of the two components, we suggest that 
IC 443 has a plasma structure consisting of a central hot 
region surrounded by an envelope of a cool and denser 
plasma, that is, the typical structure of a SNR predicted 
by the Sedov- Tayler solution. In addition, we have de- 
tected strong line emission from H-like Ne, Si and S. The 
line intensity ratio of H-like Ka to He- like Ka of S requires 
a temperature of 1.5 keV, which is significantly higher than 
the continuum temperature of 1.0 keV. The same can be 
concluded for Si. Therefore, we suggest that the 1.0 keV 
plasma is "overionized." 

In order to produce an overionized thermal plasma, 
— ™&)gas cooling rate must be higher than the recombina- 
tion rate of ions. Strong cooling due to thermal conduc- 
tion from the hot (1.0 keV) remnant center to the cooler 
(0.2 keV) exterior could cause this. 

IC 443 is the first SNR in which evidence for overi- 
onization has been found. The process of thermal con- 
duction should arise in all SNRs. Therefore we expect 
that observations of other evolved SNRs (such as mixed- 
morphology types) should detect further evidence of ove- 
rionized plasma. The apparent absence of overionized 
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gas in the similar remnant W44 suggests that the inte- 
rior plasma in other remnants may evolve differently from 
that in IC 443. 

We thank Yasunobu Uchiyama for his thoughtful com- 
ments and suggestions. The Digitized Sky Survey was pro- 
duced at the Space Telescope Science Institute under U.S. 



Government grant NAG W-2166. The images of these sur- 
veys are based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain and the 
UK Schmidt Telescope. The plates were processed into 
the present compressed form with the permission of these 
institutions. 
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R.A. offset (aromln) 



Fig. 1. — Fields of view of ASCA GIS and SIS used in the current analysis of IC 443 (20' radius region for GIS and whole SIS field of view) 
superposed on the Palomar Observatory Sky Survey red image. Offsets indicated are in arcminutes relative to the adopted center at R.A.= 
06 h 17 m 20 s , Dec.=22°40'50" (J2000). 
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R.A. offset (arcmln) R.A. offset (arcmln) R.A. offset (arcmln) 

Fig. 2. — (a): Vignetting- and exposure-corrected GIS 0.7-10.0 keV band image with both SIS and GIS FOV boundaries. GIS softness 
ratio (-Fb.7-1.5 kev/^1.5— 10.0 kov) contours over (b): GIS 0.7-10.0 keV band image, and (c): Palomar Digital Sky Survey image. The scale 
of softness-ratio contours is linear and their levels are 60%, 70%, 80%, 90% of the peak value. Grayscales of the GIS and sky survey images 
are linear and log respectively. Offset center of each figure is R.A.=06 h 17 m 20 s , Dec.=22°40'52". The contours arc limited at the FOV edge 
of the GIS along the northwest boundary, indicating the remnant extends further. The spot-like structure seen at (2,-20) in the GIS image is 
the non-thermal component so-called HXF (Kcohanc ct al. 1997; Bocchino & Bykov 2000; Olbcrt et al. 2001). 




-75 -10 -5 5 10 15 -15 -10 -5 5 10 15 

R.A. offset (arcmin) R.A. offset (arcmin) 



Fig. 3. — Exposure- and vignetting-corrected SIS images of (a): 0.5-1.0 keV, and (b): 1.0-2.0 keV bands. Two boxes are spectra extraction 
regions named "North" and "Center". Offset center of each figure is R.A.=06 h 17 m 20 s , Dec.=22°40'52". 
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Fig. 4. — SIS spectra of Center and North regions, (a): The left panel shows the Center and North spectra in black and gray respectively. 
The spectrum of the North region is re-normalized to ten times the original, (b): The right panel shows the intensity ratio of Center to North 
without any re-normalization. Below 1.4 keV, the ratio decreases proportionally with energy. 
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Fig. 5. — Ratios of the fitted Gaussian fluxes of H-like Ka to that of He-like Ka for Mg, Si, and S. Black squares and circles are those 
data from the North and Center respectively. This ratio represents the ionization degree of each element. Dashed line represents the ratios 
of H-likc Ka to that of He-like Ka with each ionization temperature calculated based on Mewe, Gronenschild, & van den Oord (1985). All 
data are on different lines of ionization temperature. 




Energy (keV) Energy (keV) 

Fig. 6. — SIS spectra of Center (left panel) and North (right panel) regions with best-fit two plasma and Gaussian models. Dashed, 
dash-dotted, and dotted lines represent the 0.2 keV plasma, 1.0 keV plasma, and additional Gaussians respectively. The lower panels show 
the residuals of the fit. 
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Fig. 7. — SIS spectra of Center (left panel) and North (right panel) regions with the best-fit model of Raymond-Smith (dashed line) and 
three narrow Gaussian components (dotted lines) with 7Vjj=7.4x 10 21 cm~ 2 at energies of 2.2-6.0 keV. The lower panels show the residuals 
of the fit and the minimum values of reduced \ 2 are 85/70 and 53/59 in Center and North respectively. 




Fig. 8. — Confidence contourT of [tk<eKi|iization temperature (T z ) to the continuum temperaturT^T|jc^\t|ie SIS spectra in Center (left 
panel) and North (right panel). fhe L confidence levels are 99%, 90%, and 67%. The contours arc above the dashed line which implies that 
the ionization temperature is higher than the continuum temperature in each figure. 
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Table 1 

asca observation logs of ic 443 







Observation Phase 




Parameter 


PV 


AO-1 


AU-D 


Sequence Number 


10011010, 10011020 51023000 


56057000 


R.A. (J2000) 


06 h 17 m 26 s 


06 h 17 m 18 s 


06 h 17 m 34 s 


Dec. 


22°40'48" 


22°26'45" 


22°20'22" 


Date 


1993 Apr 14-15 


1994 Mar 9-10 


1998 Mar 25-27 


observation mode 








SIS 


4CCD 


4CCD 


2CCD 


CIS 


standard-PH 


standard-PH 


standard-PH 


effective exposure-time 








SlSO/SlSl 


35.4/36.6 ksec 






GIS2/GIS3 


39.7/39.7 ksec 


35.4/35.4 ksec 


37.5/37.5 ksec 



Table 2 

Best-fit parameters of the two plasma model 



North Region Center Region 



model 


parameter(unit) 


value 


(90% c.r.) 


value 


(90% c.r.) 


wabs 


A H (10 21 cm- 2 ) 


7.3 


(6.7-7.9) 


7.5 


(6.9-8.4) 




/cT(kcV) 


0.17 


(0.15-0.19) 


0.18 


(0.16-0.21) 


vgnei 


n c t(10 n cm" 3 s) 


1.0 


(0.7-1.5) 


1.0 


(0.7-1.5) 


(low T) 


< kT >(keV) 


0.44 


(0.39-0.53) 


0.51 


(0.46-0.59) 




constant factor a 


185 


(117-287) 


48 


(24-87) 


vraymond 


fcT(kcV) 


1.09 


(1.01-1.19) 


1.06 


(1.01-1.11) 


(high T) 


emission measure 13 


0.038 


(0.030-0.046) 


0.056 


(0.048-0.064) 




O 


0.03 


(0.02-0.04) 


0.07 


(0.05-0.10) 




Nc 


0.09 


(0.08-0.11) 


0.12 


(0.10-0.16) 




Mg 


0.11 


(0.09-0.14) 


0.25 


(0.19-0.35) 


Abundance 


Si 


0.42 


(0.34-0.52) 


0.38 


(0.32-0.47) 




S 


0.43 


(0.35-0.54) 


0.43 


(0.38-0.49) 




Ar 


0.64 


(0.39-0.89) 


0.57 


(0.41-0.72) 




Fe 


0.07 


(0.05-0.09) 


0.17 


(0.14-0.20) 


NeX Ka 


observed energy (kcV) 


1.012 


(1.006-1.018) 


1.021 


(1.018-1.025) 




flux c 


21 


(18-26) 


29 


(23-39) 


NeX Kf3 


observed energy (keV) 






1.223 


(1.215-1.231) 




flux c 






3.1 


(2.3-4.1) 


Si XIV Ka 


observed energy (keV) 


2.015 


(2.002-2.011) 


2.011 


(2.002-2.018) 




flux c 


1.1 


(0.8-1.4) 


1.3 


(1.1-1.6) 


S XVI Ka 


observed energy (keV) 


2.664 


(2.641-2.679) 


2.649 


(2.630-2.665) 




flux c 


0.31 


(0.20-0.42) 


0.53 


(0.43-0.64) 




X Vd.o.f 


160/108 




226/116 





a This constant factor present the ratio of the VGNEI emission measure to the VRAYMOND one in XSPEC. 
b The unit is And 2 x 10 14 cm -5 where d is a distance to IC 443. 
c Thc unit is 10~ 4 photons cm~ 2 s _1 . 

Note. — Each abundance is relative to the solar value and those of two plasma components linked together. 
The 90% confidence range (c.r.) for temperature, ionization timescale, abundances, and line parameters are given in 
brackets. 

The line width of each Gaussian profile is fixed to 0. 



